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We study the Topological Casimir effect, in which extra vacuum energy emerges as a result of 
the topological features of the theory, rather than due to the conventional fluctuations of the phys- 
ical propagating degrees of freedom. We compute the corresponding topological term in quantum 
Maxwell theory defined on a compact manifold. Numerically, the topological effect is much smaller 
than the conventional Casimir effect. However, we argue that the Topological Casimir Effect is 
highly sensitive to an external magnetic field, which may help to discriminate it from the conven- 
tional Casimir effect. It is quite amazing that the external magnetic field plays the role of the 6 
state, similar to a 9 vacuum in QCD, or 9 — n in topological insulators. 



I. INTRODUCTION. MOTIVATION. 

The nature of the conventional Casimir energy is well- 
understood by now: the effect is due to the vacuum fluc- 
tuations of physical photons, which have slightly different 
propagating properties in the presence of boundaries in 
comparison with infinite Minkowski space. Essentially, 
the electromagnetic modes get modified as a result of 
nontrivial boundary conditions (BC). This tiny devia- 
tion leads to the well known expression for the Casimir 
energy 
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where a is the separation distance between the two con- 
ducting plates with size L. This extra energy gives rise 
to an attractive force per unit area (vacuum pressure) [1] 
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Today the Casimir force (2) has been measured [2], con- 
firming Casimir's basic idea [1]. 

Most of the research devoted to the computations of 
the Casimir effect deals in fact with the study of a simple 
scalar field. In this case, one can use a number of dif- 
ferent technical tools to proceed with computations [3]. 
Furthermore, it is commonly accepted that the contribu- 
tion to energy and pressure due to the physical Maxwell 
photons can be obtained directly from expressions for 
a massless scalar field by multiplying the correspond- 
ing scalar expressions by a factor of two, which accounts 
for the number of independent polarization states. This 
commonly accepted receipt follows from analysis of the 
conventional quantization of Maxwell theory in Coulomb 
gauge when unphysical degrees of freedom are completely 
removed from the system and only physical transverse 
polarizations remain in the description. This is exactly 
the procedure which leads to the well-known equations 
quoted above, (1) and (2). 

We will explicitly show, however, that this commonly 
accepted receipt is not quite correct in general when the 
theory is formulated on a compact manifold. In this case, 
there are extra contributions to the vacuum energy re- 



sulting from a new type of fluctuations due to nontriv- 
ial topological features of the gauge fields, which do not 
have counterparts in scalar field theory. These new con- 
tributions to the energy E top do depend on geometrical 
and topological properties of the system, such that the 
deviation E TC = (-Etop - -^Minkowski) 7^ leads to ad- 
ditional observable forces, similar to (2). The study of 
these new topological contributions is the main subject 
of the present work. 

A simple way to get some feeling on the nature of these 
new topological contributions is to study Maxwell the- 
ory in two dimensions, which is essentially the Schwinger 
model without fermions. As is well known, the Maxwell 
theory in two dimensions is empty, since there is no phys- 
ical propagating degree of freedom. Still, there are non- 
trivial topological sectors in the model which eventually 
lead to the emergence of the so-called 6 vacuum state. 
The construction of these topological sectors is sensitive 
to the size of the system. Therefore, it is not a surprise 
that the partition function will be also sensitive to the 
system's size, including finite size along the compactified 
time direction j3, corresponding to a finite temperature 
T = 1/(9. We will elaborate on this example with great 
details in section II using the Hamiltonian as well as Eu- 
clidean path integral approaches to explain the nature 
of topological vacuum fluctuations. We also elaborate on 
the physical "reality" of these vacuum fluctuations in Ap- 
pendix A where we compute some important observables, 
such as topological susceptibility and entropy, which are 
entirely saturated by purely topological vacuum fluctua- 
tions. 

We also note here that the extra contributions to the 
Casimir energy on a compact manifold have been dis- 
cussed previously in literature, see recent work [4] and 
references on earlier papers therein. Our goal here is 
to discuss some key elements of these new topological 
terms in a physical and intuitive way, rather than in a 
formal mathematical framework. Furthermore, we will 
discuss the relation between the 9 states and the physi- 
cal realization of these states by placing the system into 
a uniform external magnetic field. We speculate that a 
high sensitivity of the extra terms to the applied external 
magnetic field might be a key element which could allow 
one to measure these novel types of vacuum fluctuations 
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in real experiments. 

To conclude this Introduction, we want to make a few 
comments on the title of this work and the term "Topo- 
logical Casimir Effect" (TCE) which will be frequently 
used in the text below. In the given context, this term 
was introduced in [5] to emphasize that a new contribu- 
tion to the energy of the system may emerge even when 
not a single propagating degree of freedom is present in 
the system, as we mentioned above. These contributions 
occur as a result of the modification of topological fluc- 
tuations in the presence of the boundaries. They origi- 
nated from non-dispersive (contact) terms in the correla- 
tion functions and cannot be expressed in terms of phys- 
ical propagating degrees of freedom, in contrast with the 
conventional Casimir effect (CE) where physical massless 
photons are responsible for behaviour (1). The TCE rep- 
resents the portion of the energy which does not have its 
counterpart in massless scalar field theory. 

The structure of our presentation is as follows. In the 
next section, we review the relevant parts of the two di- 
mensional Maxwell "empty" theory which does not have 
any physical propagating degrees of freedom, but does 
show nontrivial topological features. We study this sys- 
tem using the Euclidean path integral approach as well 
as the Hamiltonian formalism. In section III we gener- 
alize our construction to four dimensional Maxwell the- 
ory. Numerical estimates in this case suggest that TCE is 
generally much smaller than conventional CE in normal 
circumstances. However, in section IV we advocate an 
idea that the effect is highly sensitive to a weak uniform 
external magnetic field. It is very similar to a construc- 
tion of the so-called 9 states in QCD. Finally, in section V 
which is our conclusion, we comment on some profound 
consequences the Topological Casimir Effect may have 
for cosmology. Furthermore, we advocate the idea that 
an experimental study of the Topological Casimir Effect 
in a laboratory might be considered as an investigation 
of the most intricate properties of the cosmological vac- 
uum dark energy observed in our universe. In Appendix 
A we argue, using an "empty" two dimensional Maxwell 
model, that the topological vacuum fluctuations are very 
real and very physical and must be taken into consid- 
eration to satisfy some important consistency conditions 
such as the Ward Identities. 



II. MAXWELL THEORY IN TWO DIMENSIONS 

The 2d Maxwell model has been solved numerous times 
using very different techniques, see e.g. [6-8] for a review. 
We have nothing new to say here. Our goal is in fact 
quite different: we want to review this "empty" model 
by emphasizing some elements which will be crucial for 
our discussions of the Topological Casimir Effect in four 
dimensions. 



A. Hamiltonian framework 

We consider 2d Maxwell theory defined on the Eu- 
clidean torus iS* 1 x S 1 with lengths L and /3 respectively. 
In the Hamiltonian framework we choose a Aq = gauge 
along with diAi = 0. This implies that Ai(t) is the only 
dynamical variable of the system with E = A\. The 
Hamiltonian density, the Gauss law and the commuta- 
tion relations are 

H = l -E\ c^lphys.) = 0, (3) 

[A 1 (x),E(y)]=ih6(x-y), 

where |phys.) is the physical subspace. The Gauss law is 
satisfied only for the ^-independent (zero) mode. There- 
fore, the problem is reduced to the quantum mechanical 
(QM) problem of a single zero mode living on a circle of 
circumference L. In other words, the configurations 
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are gauge equivalent and must be identified. The fact 
that 2d Maxwell theory does not describe any physical 
propagating degrees of freedom is well known- it simply 
follows from the observation that the polarization of a 
photon must be perpendicular to its momentum. How- 
ever, such a polarization can not live in the physical space 
as there is only one spatial dimension x, which is reserved 
for momentum. The presence of a single x-independent 
mode and the absence of all other x-dependent modes 
are manifestations of the "emptiness" of this theory. 

The loop integral e J dxA\ = eA\L plays the role of 
phase 4> in the conventional QM problem for a particle 
on a circle with periodic boundary conditions. The com- 
mutation relation (3) then implies that the electric field 
E is a constant in space and that it is quantized: 

E = en neZ. (5) 

The hamiltonian H = T-LL and the corresponding eigen- 
values E n for this system are well known and are given 
by 

1 d 2 1 

ff = -_. " E n = -n 2 e 2 L. (6) 
2L dAj' " 2 W 

Consequently, the partition function for this system is 
Z(/3,i) = ^ e -^=^e-^ 2 . (7) 

The construction of the so-called 9 states is also well 
known for this system [6]. The spectrum in this case 

is shifted as follows E n (9) = A (n + tp) e 2 L, such that 
the corresponding partition function now takes the form 

Z{V 1 9) = Y j e-^ n+ ^\ (8) 

where V = (3L is the two-volume of the system. Before 
we discuss the physical meaning of the obtained results 
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in the context of our present work, we want to repro- 
duce the same partition function for the same 2d Maxwell 
theory using the path integral approach 1 . In this case, 
the interpretation of eq. (8) will be quite obvious and 
straightforward. Furthermore, it can easily be general- 
ized to four dimensional Maxwell theory defined on a 
compact manifold. 



B. Euclidean Path Integral Approach 

We follow [8] and introduce the classical "instantons" 
in order to describe the different topological sectors of the 
theory which are classified by integer k. The transitions 
between different topological fc-sectors are described by 
these "instantons" . The corresponding configurations on 
the two dimensional Euclidean torus with total area V = 
L/3 can be described as follows [8]: 
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(9) 



where Q = w-E is the topological charge density and 



J d 2 x Q(x) = ^- J d 2 x E(x) = k 



(10) 



is the integer- valued topological charge in the 2d U(l) 
gauge theory, E(x) = OqAi — diA is the field strength 2 . 
The action of this classical configuration is 
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(11) 



This configuration corresponds to the topological charge 
k as defined by (10). The next step is to compute the 
partition function defined as follows 

2(6) = [vA^e~^ d2xE2+l £f d2xE . (12) 

All integrals in this partition function are gaussian and 
can be easily evaluated using the technique developed in 
[8]. The result is 
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where the expression in the exponent represents the clas- 
sical instanton configurations with action (11) and topo- 
logical charge (10), while the factor in front is due to 



the fluctuations. The computation of this pre-exponent 
factor is reduced to a conventional quantum mechani- 
cal (QM) problem as the fluctuating field is in fact x- 
independent in the Aq = gauge, as mentioned in sec- 
tion II A. Therefore, the expression for the pre-exponent 
is 
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(14) 



A simple way to evaluate this path integral is to rescale 
the A\ field according to its natural dimensionality A\ = 
ai where the dimensionless variable < a\ < 1 fluc- 
tuates inside a unit interval according to (4). In terms of 
this rescaled field problem, (14) is reduced to a standard 



expression for a free particle with mass m 
that 



m 
2^/3 
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which is precisely the pre-exponent factor in formula (13). 

While expressions (8) and (13) look differently, they 
are actually identically the same, as the Poisson summa- 
tion formula states: 
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(16) 



see [9] with detailed discussions on the relation between 
Hamiltonian formalism and the path integral approach. 



C. Interpretation 

The crucial observation for our present study is that 
this naively "empty" theory which has no physical prop- 
agating degrees of freedom, nevertheless shows some very 
nontrivial features of the ground state related to the 
topological properties of the theory. These properties are 
inherent features of the gauge theories and do not have 
counterparts in conventional scalar field theories. Rather 
these new properties are related to the presence of differ- 
ent topological sectors in the system, which we refer to as 
the "degeneracy" of the ground state, for short 3 . We in- 
terpret the nontrivial properties of the partition function 
(12,13) in this "empty" model as a result of tunnelling 
between these "degenerate" winding \n) states. These 



1 Here we use a Euclidean metric where path integral computations 
are normally performed. 

2 One should not confuse the electric field from the Euclidean for- 
mulation (9) with an E field (5) computed in Minkowski space. 
In the former case it is an unphysical complex configuration sat- 
urating the path integral, while in the latter it is a "real" phys- 
ical fluctuating electric field in the Hamiltonian formalism in 
Minkowski space. 



3 Not to be confused with the conventional term "degeneracy", 
when two or more physically distinct states are present in the 
system. In the context of this paper the "degeneracy" im- 
plies the existence of winding states \n) constructed as follows: 
T\n) = |n + 1). In this formula the operator T is the large gauge 
transformation operator which commutes with the Hamiltonian 
[7~, H] = 0. The physical vacuum state is unique and constructed 
as a superposition of \n) states as follows \9) = exp(in#)|n). 
In the path integral approach, the presence of n different sectors 
in the system is reflected by summation over k S Z in eq. (12,13). 
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tunnelling processes are happening all the time and the 
intensity of tunnelling is determined by the topological 
charge (10) and the size V of the compact manifold (11). 
A typical value of the topological charge k which sat- 
urates the series (13) in the large volume limit is very 
large, k ~ Ve 2 V 3> 1. 

It is different from the conventional tunnelling in QM 
in that the tunnelling in our system corresponds to 
a transition between one and the same physical state, 
whereas that in QM describes a transition between phys- 
ically distinct states. More specifically, the tunnelling 
in our case occurs between the winding |n) states which 
are connected by large gauge transformations. Therefore, 
they correspond to one and the same physical state. 

The key for our present work is the observation that 
the properties of these tunnelling processes are sensitive 
to the size of the system. In different words, the addi- 
tional energy associated with these tunnelling processes 
is different for systems with different sizes and shapes. A 
direct manifestation of this sensitivity (when it is gener- 
alized to 4d case as we discuss below) is the emergence of 
the Topological Casimir Effect (TCE), when the vacuum 
energy and pressure depend on the size of the compact 
manifold on which the theory is defined. 

Is this extra energy physical? Our ultimate answer is 
"yes". We refer to Appendix A where we present some 
arguments suggesting that the extra energy related to 
the tunnelling processes in the "empty" theory can not 
be removed by any redefinition of observables. It must 
be present in the system for consistency of the theory. 
In particular, the Ward Identities can not be maintained 
without these tunnelling contributions, see Appendix A 
for details. Essentially, this extra contribution is pre- 
cisely the source of violation of a commonly accepted 
(but generally wrong) receipt that the Casimir effect due 
to Maxwell photons could be obtained by multiplying the 
corresponding scalar expressions by a factor of two. 



III. TOPOLOGICAL CASIMIR EFFECT IN QED 
IN FOUR DIMENSIONS 

The topological structure of the gauge field in 2d can 
be easily generalized to higher dimensions. In four space- 
time dimensions, we can devise boundary conditions that 
give rise to very similar instanton-like configurations, 
with precisely the action (11) as found in 2d Maxwell 
theory on the torus. In this section we show that these 
topological degrees of freedom are completely decoupled 
from the propagating physical photons. Furthermore, the 
corresponding quantum fluctuations do not depend on 
the properties of the topological sectors (due to the lin- 
earity of the Maxwell equations), and can be treated in 
the conventional way. As a result we are able to focus on 
the new contributions and compare them to the conven- 
tional Casimir effect from literature. We shall see that 
the topological Casimir effect is strongly suppressed on 
a Euclidean 4-torus where one of the spatial dimensions 



is much smaller than the others. In this case the well 
known formulae (1) and (2) are recovered. However, the 
main goal of this work is to study precisely those novel 
contributions which are sensitive to the system size. 

A. Decoupling of the topological and conventional 
parts 

To construct a theory defined on a Euclidean 4-torus, 
we consider a system with a box of sizes L\ x X L3 x f3 
in the respective directions. A torus is realized when 
we assume periodic boundary conditions on the physical 
fields in all directions, in which case we find a "degen- 
eracy" of the vacuum state, just like in section II: by 
making a loop in the xy-plane, the field can pick up a 
phase corresponding to a large gauge transformation. In 
fact, the current physical setup with the torus can be de- 
scribed in a manner similar to the usual scenario with the 
parallel plates (i.e. with plates lying parallel to the xy- 
plane and separated by a distance in z), but demanding 
the non-trivial BCs given above. Working in Euclidean 
space and adopting the Lorentz gauge, it is simple to find 
a 4d generalization of the instanton potential from sec- 
tion II that satisfies these boundary conditions. The 4d 
instanton potential is given by 

,,, / irk irk \ 

where k is the winding number that labels the topologi- 
cal sector, and L\, L2 are the dimensions of the plates in 
the x and y-directions respectively, which are assumed to 
be much larger than the distance between the plates L3. 
This classical configuration satisfies the periodic bound- 
ary conditions up to a large gauge transformation, and 
provides a topological magnetic flux in the z-dircction: 

B top = V x A top = f 0, 0, -^j- J , (18) 

in close analogy with the 2d case (9). The Euclidean 
action of the system becomes 

l - J d 4 x^E 2 + (B + B top ) 2 y (19) 

where the integration is over the Euclidean torus L\ x 
L2 x L3 x (3 and E and B are the dynamical quantum 
fluctuations of the gauge field. These terms were not 
present in the 2d model, but must here be taken into 
account due to the presence of real propagating physical 
photons. We find that the action can be easily split into 
the sum of a topological and a quantum part, because of 
the vanishing cross term 

J d 4 x B ■ B top = -M- J d 4 x B z = (20) 

Here the fact is used that the magnetic portion of quan- 
tum fluctuations in the z-direction, represented by B z = 



5 



d x A y — d y A x , is a periodic function because A is periodic 
over the domain of integration. As a result, there is no 
coupling between the conventional quantum fluctuations 
described by photons with physical polarizations and the 
classical instanton potential (17), (18). Furthermore, the 
quantum fluctuations due to photons are not sensitive to 
the topological sector k of the theory, and therefore they 
decouple from the classical fc-instanton contribution. Fi- 
nally, the quantum fluctuations from photons must be 
computed in a box with size LiL 2 rather than in infinite 
space along x, y directions. The corresponding correc- 
tions, in principle, can be computed. Technically the 
computations would be quite tedious as they require the 
operation with Green's functions defined on a finite man- 
ifold rather than in the infinite space. The computations 
can be performed for the trivial k — topological sector 
as the corresponding corrections are independent of k. 
These contributions are expected to produce some cor- 
rections ~ (1 + jp) to formula (2). However, we shall 
not elaborate on these terms in the present work. It is a 
part of the conventional partition function Z computed 
for trivial topological sector k = 0. 

The main lesson from the previous discussion is that 
the conventional quantum fluctuations are not sensitive 
to the topological sectors k as a result of linearity of the 
Maxwell equations. Therefore, they can be treated in 
completely separate ways, which greatly simplifies our 
analysis. For the partition function we can now write: 
Z = Z x Z top . The conventional part Z is well- 
studied for toroidal BCs at finite temperatures [3]. It 
is k- independent, so it will not be elaborated here. In 
the rest of this section we will study the behaviour of 
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Ztop- We shall see that in the limit when =^ 
the partition function related to the topological effects 
yield Z top = 1, and we recover the conventional Casimir 
effect (1), (2) which is computed from Zq. However, we 
shall see that a number of novel and unusual features will 
emerge in the system when L\ , L 2 are large but remain 
finite, which is precisely the main subject of our studies. 



B. Computing the topological pressure 

The system of parallel plates is related to 2d Maxwell 
theory by dimensional reduction: taking a slice of the 4d 
system in the xy-plane will yield precisely the topologi- 
cal features of the 2d torus. Assuming that L3 is much 
smaller than L\ and L 2 , the additional dimensions do not 
contribute toward Z top as we noted above. Instead, the 
quantum corrections slightly modify Zq as they do not 
depend on the topological sectors k, and can be factored 
out from Z top . With this set up, the classical action for 



configuration (18) takes the form 
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while the topological partition function becomes: 
27r/3L 3 
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where the 2d electric charge entering eqs. (11), (13) is 
expressed in terms of the 4d electric charge as follows 4 , 
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In this section we consider 9 = 0. More discussion on 
this matter is found in section IV. 

One should note that the dimensional reduction which 
is employed here is not the most generic one. In fact, 
one can impose a non-trivial boundary condition on every 
slice in the 4d torus. However, the main goal of this work 
is not to classify the most generic BC, but to discuss the 
physical properties for the simplest possible case (18), i.e. 
when a nontrivial BC is imposed on a single slice, while 
keeping the trivial periodic BC for other slices. 

With this objective in mind, it is useful to introduce 
the dimensionless parameter 



t = 2pL 3 /e 2 L 1 L 2 



(24) 



such that the partition function Z top can be written in a 
very simple form: 
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(25) 
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where the Poisson summation formula (16) is used again. 
Our normalization of the partition function Z top is such 
that in the limit LiL 2 = 00 the topological portion of 
the partition function Z top = 1 so that we recover the 
conventional Casimir effect (1), (2) which is encoded in 
Zq. The simplest way to check our normalization is to 
take the limit r — > using the right hand side of eq. 
(25) when a single term with n — contributes. It corre- 
sponds to very large instanton numbers k ~ t -1 / 2 — > 00 
saturating the original series (22). 

From Z top , we can calculate any thermodynamic prop- 
erty of the system, like the topological pressure between 
the plates 
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Ptop ~ pL,L 2 dL 3 lnZtop - 



(26) 



In the asymptotic limit where r <C 1 one can use the 
dual representation of Z top encoded by the Poisson re- 
summation formula (25) to find that 
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T < 1. 



(27) 



4 Note that the units are consistent. Indeed, in 4d, e 2 ~ a is the 
dimensionless fine-structure constant. However, in 1+1 dimen- 
sions the QED coupling constant has units of {length) -2 . 



In this case, the topological pressure is exponentially 
suppressed, and when compared with the conventional 
Casimir pressure (2) it is clear that the topological effect 
is too small to measure experimentally 5 . 

A few comments are in order. First of all, one can 
explicitly see that our original instanton formula (22) is 
consistent with our interpretation: that the additional 
energy and pressure due to the topological features of the 
system is the result of tunnelling events between different 
winding states, as discussed in section II C. Indeed, a non- 
analytical structure of eq. (22) with respect to coupling 
constant exp(— 1/e 2 ) represents the typical behaviour for 
a tunnelling process. It is quite fortunate that the Pois- 
son re-summation formula (25) allows us to analyze both 
regimes, at large as well as small r. Secondly, even in 
this simple t <C 1 case one can explicitly see that the 
sign of the effect is opposite to the conventional Casimir 
effect (2). This correction leads to repulsive rather than 
attractive forces. This "wrong" sign is a typical manifes- 
tation of the topological fluctuations, in contrast with the 
conventional vacuum fluctuations of photons with phys- 
ical polarizations. See some additional comments on a 
"wrong sign" in Appendix A. 

While t<1 can be examined analytically, it is more 
interesting to study a system with r ~ 1 where this effect 
could be sufficiently large. We can satisfy the conditions 
Li,L,2 S> P,Ls, in order to use the dimensional reduc- 
tion, and at the same time still achieve r ~ 1 because the 
small parameter e 2 enters the denominator in eq. (24) in 
the definition for r. In this case we have to resort to 
numerical approximations, since there is no closed form 
for the partition sum and pressure. In Figure 1, a nu- 
merical plot of Ptop is shown for this regime. There is 
a large peak around r ~ 0.4 when pressure, measured 
in units L 2^2 e 2 > has an order of 1. The relative mag- 
nitude between the maximum topological pressure and 
the conventional Casimir pressure (2) using parallel ideal 
conductors at low temperature is thus approximately 

^ |P top | ^ 480£| ^ 120 Lj 
max ~ \P\ ~ L\L\e^ ~ ifa ' L\L\ ' [ ' 

This ratio (even at its maximum at r ~ 0.4) is very small 
in a typical Casimir experiment setup with L\, L2 S> £3, 
in spite of the large numerical factor in front of formula 
(28). As we mentioned previously, the power-like correc- 
tions ~ L\lL\,L\jL\ are also expected to occur in Zq 
resulting from the conventional vacuum fluctuations of 
physical photons. We expect these conventional correc- 
tions to be even smaller as they can not contain a para- 



5 The conventional Casimir pressure in eq. (2) does not account 
for the thermal correction, and is computed for a system with 
slightly different boundary conditions (metallic instead of peri- 
odic). However, the boundary conditions only change the pres- 
sure by a constant of order unity and the thermal correction 
is negligible for low temperatures, see [3]. Also note that even 
though t<1, wc are still in the low-temperature regime. 
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metrically enhanced factor 1/e 2 that is a unique feature 
of the topological vacuum fluctuations. 



Topological Casimir Pressure 
1-1 



P 




1 



FIG. 1. The topological pressure on the 4d system of par- 
allel plates as a function of r s 2j3Lz/ e 2 L\L2- Pressure is 
measured in units 2 2 2 2 . 

To conclude this section, we find that there is a small, 
but very real, contribution to the Casimir effect that is 
purely due to topological features of the system. When 
QED is defined on a compact manifold such as a 4-torus, 
one needs to take into account the tunnelling processes 
which occur between the topologically inequivalent (but 
physically identical) winding states. These topological 
transitions are described in terms of integer magnetic 
fluxes (18). It is not surprising that the effect is expo- 
nentially small in normal circumstances (27). The effect 
remains very small (28) even at r ~ 1. Still, there is 
a hope to make it measurable by studying the Topo- 
logical Casimir Effect (TCE) in the presence of some 
external magnetic field. We shall observe a high sen- 
sitivity of TCE to applied weak external magnetic field. 
This should be contrasted with conventional Casimir ef- 
fect (1), (2) which can not be sensitive to external fields 
as vacuum photon fluctuations do not couple to exter- 
nal fields (since the Maxwell equations are linear). This 
topic is precisely the subject of the next section. 

IV. 9 VACUA AND EXTERNAL MAGNETIC 
FIELDS 

Now it is interesting to place our system into a region 
with a weak external magnetic field B% xt along the z- 
direction. The idea behind of this construction is that the 
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external magnetic field B% xt will interfere with the inte- 
ger topological flux (18) describing the tunnelling events. 
It is expected that such interference may skew the sum- 
mation over the topological sectors, similar to the ac- 
tion of the so-called 9 parameter (12), (13). As we shall 
demonstrate below, this is indeed what happens in our 
simple case considered in section III. In different words, 
we claim that by adding a constant magnetic field to 
the previous setup, an effective non-zero 6 e f / parameter 
emerges in the system. The crucial point here is that we 
introduce this parameter which can be externally var- 
ied. By studying the corresponding responses to 9 e ff 
variation, it gives us some hope that while the TCE is 
numerically very small (28), it is nevertheless very sensi- 
tive to a weak magnetic external field (in contrast with 
conventional Casimir effect (2)), and hopefully it can be 
eventually measured due to this sensitivity. 

To construct a system as such, in addition to the topo- 
logical flux through the xy-plane, we apply a real physical 
constant magnetic field B% xt — d x A ext — d y A% xt parallel 
to the z-direction (perpendicular to the xy-plane). The 
total B z field in the Euclidean metric is thus modified as 
follows: 



B, 



Q1 _|_ Q to V _|_ B ext 



(29) 



where the total field decomposition consists of the same 
instanton potentials A l ° p as in eq. (17), the external 
magnetic field potential A ext given above and the quan- 
tum fluctuations Af. around them. 

The only difference from the previous construction is 
the additional external constant magnetic field. Note 
that the quantum fluctuations still decouple from the 
classical and external fields, similar to eq. (20), due to 
the periodicity of quantum fluctuations over the domain 
of integration, 



AS = (Bt 



2irk 
LiL 2 e 



) • / d^xBl = 0. 



(30) 



The remaining part of the action is quadratic and thus 
path integration can be performed. The same calculation 
from the previous section follows and the partition func- 
tion separates into a classical portion, which describes 
TCE, and a quantum portion that corresponds to the ef- 
fect of photons in 4D, i.e. the well known Casimir effect. 
It is important that the conventional quadratic term rep- 
resenting the photon fluctuations does not depend on the 
topological sector k, nor on the external magnetic field. 
It is described exclusively by Zq, as before. 

The new element here is that the external field couples 
to the instanton potential. Therefore, the topological 
partition function now takes the form, 



Z top (T,8 ef f) = ^/ttt y^exp 



kez 



-7r 2 r k + 



°eff 

2vr 



(31) 



where we introduced the effective theta parameter 

9 eff = Bl xt LiL 2 e (32) 



proportional to the external magnetic flux through the 
xy-plane in this particular system. It is clear from the 
partition sum (31) that a non-zero effective Q e ff skews 
the summation over topological sectors similar to the 2d 
example given by (12), (13). It is also clear that 9 e ff — 
27rm corresponds to integer flux m through the xy-plane, 
which obviously can not modify the system, such that 
Ztop(T,0 e f f) is 2ir periodic in 9 e ff. 

In what follows, we also need a "dual" representation 
for Z top (T,9 e ff) which is obtained by applying the Pois- 
son re-summation formula (16) 



Z top (T,e ef f) = J2 ex P 



+ i Oeffn 



(33) 



In representation (33), it is obvious that 9 e ff being ex- 
pressed in terms of the external magnetic field (32) can 
be thought of as a fundamental 9 parameter. How- 
ever, the corresponding "instanton charge" n which nor- 
mally enters with 9 is not the same magnetic flux k from 
our original construction (18) with classical action (21). 
Rather, it is some "dual" configuration with classical ac- 
tion ~ T -1 . 

We note that Z top {T,9 e f /) is properly normalized in 
the limit L\.Li — s- oo which corresponds to r —¥ 0. In 
this case, only a single term with n — in (33) survives, 
leading to the desired normalization Z top (r — > 0, e //) — 
1. Therefore, all conventional formulae for the Casimir 
effect determined by Zq are recovered in this limit as Z 
is factorized Z = Z x Z top for our system. 

Now we are in position to calculate the topological 
Casimir pressure contribution from free energy at finite 
temperature, by inserting the partition function into eq. 
(26). Like in section III B, the topological pressure has no 
closed form. For the limit where r <C 1, we may obtain 
an asymptotic expansion. Using the dual representation 
(33) and keeping the leading order terms, we arrive at: 



Pi 



top 



cos(6» e// )exp(-l/-r), t<1 (34) 



which reduces to our previous formula (27) in the ab- 
sence of the external magnetic field. As expected, the 
oscillatory effect with respect to 9 e ff is present, but be- 
comes exponentially suppressed because the tunnelling 
amplitudes naturally diminish in this limit. 

In a more general case when r ~ 1, we have to use 
some numerical methods as asymptotic analysis is no 
longer sufficient. A numerical plot is shown in Figure 
2, such that the variation with 9 e ff is manifest. The 
pressure is clearly oscillatory with respect to 9 e f j and its 
local extrema are attained at mr where n € Z. Thus, by 
altering the magnetic flux, the topological Casimir pres- 
sure will also be modified accordingly. Additionally, at 
9 e ff — the pressure has a "wrong" sign, i.e. it is op- 
posite to conventional Casimir effect (2), as we already 
discussed after eq.(27). This sign changes as a function 
of the external magnetic field, as can be seen in the plot. 
Such a variation can be interpreted as the result of in- 
terference between the external magnetic field and the 
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FIG. 2. Topological Casimir Pressure is plotted in units 
L 2^i e 2 f° r different r. A clear 2tt periodicity is seen and local 
extrema are between odd and even integer multiples of -n. 



{Bind) = for 6 e ff € {nn : n £ Z}. This property is 
analogous to the well known property in QCD when V 
and CV invariance holds only for 9 = nn. 
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topological "instanton" fluctuations (18). This very spe- 
cific and distinct variation, hopefully, might be a useful 
feature to measure the effect in the future in spite of its 
strong numerical suppression (28) in comparison with the 
conventional Casimir effect (2). 

Our next step is to analyze the magnetic response of 
the system under influence of the external magnetic field. 
The idea behind these studies is the observation that the 
external magnetic field acts as an effective 9 e f / parameter 
as eqs. (32), (33) suggest. Therefore, the corresponding 
study of the magnetic response in our system is in fact 
very similar to an analysis of the expectation value of the 
topological density and topological susceptibility in other 
gauge models, where a 9 parameter enters the system in 
a similar way. See e.g. Appendix A for some discussion 
on this in 2d Maxwell theory. 

First, we compute the induced magnetic field defined 
as follows 

, _ 1 d\nZ top _ e d\nZ t0 p 

Wind)- p y ~ d0eff W 

= ^ £ (V + exp + °f±) 2 . 

As one can see from (35), our definition of the induced 
field accounts for the total field which includes both 
terms: the external part as well as the topological por- 
tion of the field. In the absence of the external field 
^gext _ ser i es j s antisymmetric under fc — > — k 

and {B in< i) vanishes. It is similar to the vanishing expec- 
tation value of the topological density in gauge theories 
when 9 = 0. One could anticipate this result from sym- 
metry arguments as the theory must respect V and CV 
invariance at 9 = 0. 

The expectation value of the induced magnetic field 
exhibits similar In periodicity from the partition func- 
tion and it reduces to triviality whenever the amount 
of skewing results in an antisymmetric summation, i.e. 



FIG. 3. A numerical plot of the induced magnetic field in 
units Ll 1 L2e as a function of 9 e ff- The oscillatory behaviour 
becomes more pronounced for large r. 

As before, the topological effects are exponentially sup- 
pressed at r < 1, as Ztop — ^ 1 with exponential accu- 
racy at t < 1. The effect is much more pronounced 
in the range where t ~ 1, see Figure 3, where we plot 
the induced magnetic field in units (L\Lie)~ x as a func- 
tion of 9 e ff. One should also remark here that the in- 
duced magnetic field defined as (35) can be thought of 
as the magnetization of the system per unit volume, i.e. 
(M) = —(Bi n d), as the definition for (M) is identical to 
(35) up to a minus sign. 

Now we turn our attention to the magnetic suscepti- 
bility, which is similar to the topological susceptibility 
reviewed in Appendix A for 2d QED. This object is V 
and CV even and does not vanish at zero external field. 
The magnetic susceptibility measures the response of free 
energy to the introduction of a source term, which is rep- 
resented in our case by B ext ~ 9 e ff. To be more precise, 
we define Xmag m a w &y which is similar to the topolog- 
ical susceptibility in Appendix A for 2d QED, 

Xm ag = J d 4 x(B z (x),B z (0)) = -± p - 2 ^P ,(36) 

where the integration is over the Euclidean torus L\ x 
L 2 x L 3 x (3. With this definition, Xmag is a dimension- 
less parameter, in contrast with 2d QED where Xeum 
has dimension (mass) 2 , and in 4d QCD where xqcd 
has dimension (mass) A . This is due to the fact that the 
topological density operator has dimension (mass) 2 in 2d 
QED and (mass) 4 in 4d QCD while our topological in- 
stanton (18) expressed in terms of the magnetic field has 
dimension (mass) 2 . Nevertheless, we opted to keep defi- 
nition (36) without inserting any additional dimensional 
parameters such as LiL 2 into formula (36), to main- 
tain the conventional definition in statistical mechanics 
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where Xmag is a dimensionless parameter (in units where 

h = c = i). 

We can represent (36) in terms of dimensionless vari- 
ables as follows, 



Xn 



2d 2 \nZ top (T,9 eff ) 



(37) 



In the limit when t«1, one can use analytical expres- 
sion (33) to conclude that Xmag ~ exp(— ^) is strongly 
suppressed. It is consistent with our expectations that 
there should not be any magnetic correlations in the con- 
ventional Casimir experimental setups. However, for r 
near the order of 1, the behaviour is quite nontrivial as 
shown in Figure 4. Note that in this case, the susceptibil- 



Magnetic Suceptibility 




FIG. 4. Numerical plot of magnetic susceptibility Xmag as 
a function of 9 e ff for different values of r. It oscillates with 
O e ff as it should, and it does not vanish even at zero external 
magnetic field at d e ff = due to topological fluctuations, 
similar to well studied cases of 2d QED and 4d QCD. 

ity is highly dependent on the external field and changes 
signs, which is extremely unusual for conventional sys- 
tems. More specifically, Xmag behaves like a weak dia- 
magnetic for 9 e ff s» nn, n even; and behaves as strongly 
paramagnetic for n odd. 

We note that this topological effect is quite distinct 
from the behaviour of the conventional Casimir effect. In 
the conventional quantization of electromagnetic fields in 
Minkowski space, there is no direct connection between 
the Casimir pressure and an external magnetic field. Any 
such coupling could only be attained through fermionic 
interactions in higher order diagrams, and it should be 
highly suppressed ~ a 2 B e 2 It /m* as the non- linear Euler- 
Heisenberg Effective Lagrangian would suggest. In a pure 
Maxwell theory, as we consider here, the photon-photon 
interactions should be trivial. However, in TCE, we see 
that the external magnetic field does have a non-trivial 
effect on physical quantities through its interactions with 
topological "instanton" fluctuations (18). 

Thus, the periodicity in all of the physical quantities 
with respect to the external magnetic field is a unique fea- 
ture of TCE. It is not found in any of the typical Casimir 



results and can serve as a clear indicator to distinguish a 
topological effect from conventional Casimir effects. 



V. CONCLUSION AND FUTURE DIRECTIONS 

We start our conclusion with the following question, 
which the readers have likely been asking themselves 
many times while reading the text. What went wrong 
with the standard and generically accepted arguments, 
which claim that the Casimir effect in Maxwell theory is 
simply twice the corresponding expression derived using 
massless scalar field theory? It is known why this ar- 
gument is not quite correct when Maxwell theory is for- 
mulated on a non-simply connected, compact manifold. 
The standard "naive" argument neglects the topologi- 
cal sectors, which are indeed absent when the theory is 
formulated on the topologically trivial Minkowski space. 
However, these sectors are present when the theory is 
formulated on a nontrivial manifold. Physics related to 
pure gauge configurations describing the topological sec- 
tors does not go away when one removes all unphysical 
degrees of freedom as a result of gauge fixing; instead, 
this physics reappears in a much more complicated form 
where the so-called Gribov ambiguities 6 emerge [10]. See 
recent paper [4] and also some previous relevant discus- 
sions [11-13]. We opted to keep some gauge freedom in 
our analysis to study these topological sectors explicitly. 

Now we can formulate the main results of the present 
paper. First of all, the physics behind of the Topologi- 
cal Casimir Effect (TCE) is quite simple: there are tun- 
nelling events in a theory formulated on a small compact 
manifold when the temperature is small but not identi- 
cally zero. These transitions can be completely ignored 
for relatively large systems, but in general, these topo- 
logical transitions interpolating between different wind- 
ing states (which correspond to one and the same unique 
physical state) do occur. The amplitudes for these transi- 
tions depend on size and shape of the system. Therefore, 
it is not really a surprise that the vacuum energy associ- 
ated with these tunnelling events depends on the size of 
the system, which ultimately implies an extra contribu- 
tion to the Casimir pressure. 

In general, the effect is numerically much smaller than 
the conventional Casimir effect with the ratio given by 
eq.(28). However, we argued that the effect is highly 
sensitive to small external magnetic fields which can serve 
as a clear indicator to distinguish TCE from conventional 
Casimir effects. 



6 A short historical remark is warranted here. The Gribov ambi- 
guities [10] were originally discussed for non-abelian gauge theo- 
ries in Minkowski space when one tries to completely remove all 
unphysical degrees of freedom in the Coulomb gauge. The corre- 
sponding Gribov ambiguities lead to strong infrared singularities 
in the consequent analysis. Another option is to deal explicitly 
with some unphysical degrees of freedom that effectively describe 
these topological sectors. 
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Our last comment is as follows. The TCE as we al- 
ready mentioned in the introduction is a very generic 
phenomenon in gauge theories. It shows the algebraic 
sensitivity to the size of the system even when the theory 
has a mass gap. See [5, 14] and many references therein 
where TCE has been tested in various models, including 
QCD lattice computations. Our comment here is that 
the observed Dark Energy (DE) in the universe might be 
a direct manifestation of the TCE as argued in [5, 14] 
and references therein. The idea is based on two key ele- 
ments. Firstly, the additional energy in Maxwell theory 
(defined on a compact manifold and discussed in this pa- 
per) is based on nontrivial topological properties formally 
expressed by the first homotopy group 7Ti[/7(l)] ~ Z. 
In four dimensions a similar structure emerges for non- 
abelian QCD where the third homotopy group is non- 
trivial 7T3[5f (3)] ~ Z. In this case one can argue that 
the system is algebraically sensitive to very large dis- 
tances in spite of the fact that the theory has a mass 
gap. The second key element is based on the paradigm 
that the relevant definition of the energy which enters the 
Einstein equations is the difference AE = (E — i?Mink)j 
similar to the Casimir effect (1), rather than the energy E 
itself. In this case the difference between the two metrics 
(expanding universe with Hubble expansion rate H and 
Minkowski space-time) as a result of TCE would lead to 
an estimate [5, 14] 

AE ~ L- X A| CD ~ (10- 3 e^) 4 , (38) 

where L is the visible size of the universe estimated as 
L~ l ~ H ~ 10 _33 eV. Estimation (38) is amazingly close 
to the observed DE value today. In fact, a comprehensive 
phenomenological analysis of this model (the so-called 
"ghost dark energy" model) has been recently studied in 
a number of papers where comparisons have been made 
with the current observational data. (See references on 
observational papers in [5, 14].) The conclusion was that 
the model (38) is consistent with all presently available 
data, and we refer the reader to the original papers on 
analysis of the observational data. 

Our comment relevant for the present study is that 
some very fascinating topological properties of the quan- 
tum vacuum may be, in principle, studied in a laboratory 
if the TCE in Maxwell theory, which is the main subject 
of the present work, can be experimentally measured. We 
conclude on this optimistic note 7 . 
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Appendix A: Why topological vacuum fluctuations 
must be real and physically observable. 

The main subject of the present work is zero point 
(vacuum) fluctuations. There has always hung a shadow 
over this question as there have always been suspicions 
that those vacuum fluctuations are not really zero point 
fluctuations, but rather can be attributed to some other 
physics. See in particular the relatively recent paper [16] 
where it has been argued that the conventional Casimir 
effect can be computed without even mentioning such a 
notion as the "vacuum" . 

The main goal of this appendix is to argue, using a 
simple exactly solvable 2d model, that the topological 
vacuum fluctuations are very real and very physical. In 
other words, we want to present a few arguments sug- 
gesting that finite contributions resulting from topolog- 
ical features of the system cannot be removed by any 
means such as subtraction or redefinition of observables. 

We start our study with the topological susceptibility 
X defined as follows, 



X 



—t lira 

47T J fc->0 



d 2 x e lkx (TE(x)E(0)) , (Al) 



where Q = ^-E is the topological charge density and 

J d 2 x Q(x) = ~ J d 2 x E{x) = k (A2) 

is the integer valued topological charge in the 2d U(l) 
gauge theory, E(x) = d\A 2 — d^Ax is the field strength. 
The x measures response of the free energy to the intro- 
duction of a source term 



La = ii 



2tt 



/ A EM 



(A3) 



The corresponding computations can be easily carried 
out as the partition function 2(6) is known exactly, see 
section II. When differentiating the partition function 
twice with respect to 6, we get a finite contribution in 
the infinite volume limit, V = f3L — > oo, i.e 
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7 In fact, the idea to test some intriguing vacuum properties rele- 
vant for cosmology in a laboratory is not a very new idea. It has 
been advocated by G. Volovik for years, see recent review [15] 
and references therein. 
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(A4) 



Is contribution (A4) physical? This question immediately 
arises because we are dealing with "empty" Maxwell the- 
ory in two dimensions, where there are no physical propa- 
gating degrees of freedom in the system. Can we redefine 
the theory to remove all such terms from consideration 
once and for all? For example, one can use a prescrip- 
tion [17] which ignores the topological sectors and leads 
to a trivial partition function Z = 1, see eqs (4.1), (4.3) 
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in [17]. Such a prescription would obviously be consis- 
tent with the conventional procedure which relates the 
Casimir effect for 4d Maxwell theory and 4d massless 
scalar field theory up to factor 2. However, such a pre- 
scription would not produce the contact term (A4) which 
must be present in the system for its consistency, as we 
shall argue below. 

The question addressed above on "physical reality" of 
(A4) is not a purely academic question. If eq. (A4) 
is treated as a physical contribution, then the partition 
function Z(/3,L,9) which leads to (A4) is also physical. 
Therefore, the same partition function Z(f3, L, 9) will also 
lead to an extra Casimir force P ~ d In Z jdL which "in 
principle" is an observable quantity. In different words, 
if (A4) is physical, then there is an extra term in the 
Casimir energy which is not related to any asymptotic 
propagating degrees of freedom and is "in principle" ob- 
servable. We present a few arguments below to advocate 
that (A4) is indeed physical. 

Our argument goes as follows. We add a massless 
fcrmion field tp to the system to arrive at the well known 
2d Schwinger model. The expression for the topological 
susceptibility in the 2d Schwinger QED model is known 
exactly [8, 14] 



Xqed 



4tt 2 



d 2 x 



d 2 (x) 



2ir 2 



K (n\x\) 



, (A5) 



where /i 2 = e 2 /ir is the mass of the single physical state in 
this model, and K (fj,\x\) is the modified Bessel function 
of order 0, which is the Green's function of this massive 
particle. The crucial observation here is as follows: any 
physical state contributes to xqed with negative sign 



Xdi 



lim > 



: E\n)(n\£E\0) 



< 0. (A6) 



In particular, the term proportional to —Ko{p-\x\) with 
the negative sign in eq. (A5) is the result of the only 
physical field of mass [i. However, there is also a contact 
term -f^ in eqs. (A5), (A4) which contributes to the 
topological susceptibility x with the opposite sign, and 
which can not be identified according to (A6) with any 
contribution from any asymptotic state. 

The first term in this formula (A5) can be easily 
recognized as the expression for xeum for 2d Maxwell 
theory (A4) which is originated from the topological sec- 
tors, and not related to propagating degrees of freedom. 

This term has a fundamentally different, non- 
dispersive nature. In fact it is ultimately related to dif- 
ferent topological sectors discussed in section II. This 
contact term must be present in the expression (A5) 
to satisfy the Ward Identity (WI) which states that 
XQED( m = 0) = 0, see [14] for the details. Without 
this contribution, it would be impossible to satisfy the 
Ward Identity because the physical propagating degrees 
of freedom can only contribute with sign (— ) to the cor- 
relation function as eq. (A6) suggests, while WI requires 
X = in the chiral limit m = 0. One can explicitly 



check that WI is indeed automatically satisfied only as 
a result of exact cancellation between conventional dis- 
persive term with sign (— ) and non-dispersive term (A4) 
with sign (+), 
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K (n\x\) 



(A7) 
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47T 2 



[1 - 1] = 0. 



The lesson we learn from this simple exercise is that the 
contact term (A4) which is saturated by the topological 
sectors is physical, and it must be present in the system 
for its consistency. 

The same contact term (A4, A5) can be also computed 
using the auxiliary ghost fields, the so-called Kogut- 
Susskind (KS) ghost, as has been originally done in ref. 
[18], see [14] for relevant discussions in the present con- 
text. This auxiliary ghost field effectively takes into ac- 
count the presence of topological sectors which lead to 
(A4). The crucial element accounting for different topo- 
logical sectors of the underlying theory, does not go away 
in KS-description. Rather, this information is now coded 
in terms of the unphysical ghost scalar field which pro- 
vides the required "wrong" sign for contact term (A4, 
A5). The contact term in this framework is precisely 
represented by the ghost contribution replacing the stan- 
dard procedure of summation over different topological 
sectors. At the same time, this unphysical ghost scalar 
field does not violate unitarity or any other important 
properties of the theory as consequence of Gupta-Bleuler- 
like condition on the physical Hilbert space, see [14] for 
the details in the given context 8 . 

Our second argument that the topological sectors must 
be taken into consideration is based on analysis of the 
entropy in the same 2d "empty" Maxwell theory. Before 
we formulate our argument, we want to make a short 
historical detour on the entropy studies in this "empty" 
model. 

It has been claimed [19] that for spins zero and one-half 
fields, the one loop correction to the black hole entropy 
is equal to the entropy of entanglement, while for a spin 
one Maxwell field, the entropy has an extra term de- 
scribing the contact interaction with the horizon. While 



8 It is important to emphasize that the KS ghost should not be 
confused with the conventional Fadeev-Popov ghost which is nor- 
mally introduced into the theory to cancel out unphysical polar- 
izations of the gauge fields. Instead, the KS ghost is introduced 
to account for the existence of topological sectors in the theory. 
A similar construction is also known for four dimensional non- 
abelian gauge theory where the corresponding color singlet field 
is called the Veneziano ghost, see [14] for references and details. 
In the four dimensional case the Veneziano ghost can not be con- 
fused with a Fadeev-Popov ghost as the Veneziano ghost being 
a singlet does not carry a color index, in contrast with Fadeev- 
Popov ghosts. The sole purpose of the Veneziano ghost is to 
saturate the contact term with the "wrong sign" in the topolog- 
ical susceptibility, similar to eq. (A4), (A5) in the 2d Schwinger 
model. 
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the entropy is a positively defined entity, the Kabat con- 
tact term is negative [19]. Furthermore, this term be- 
ing a total divergence can be represented as a surface 
term determined by the behaviour of the theory at arbi- 
trarily large distances, i.e. it obviously has an infrared 
(IR) origin. More recently, it has been conjectured [14] 
that the Kabat contact term is originated from the same 
topological gauge sectors which saturate the topological 
susceptibility (A4). Indeed, both terms have "wrong" 
signs in comparison with what physical propagating de- 
grees of freedom would produce, and both terms can be 
represented by surface integrals, see [14] for the details. 
Next step in this development was the computation of 
the entropy for the 2d Maxwell system, defined on a fi- 
nite dimensional compact manifold with size V = f3L, 
such that the IR physics can be properly treated [20], 
see also [21] with related discussions. In this case the 
expression for the entropy can be easily computed from 
the partition function (8), (13), (16) discussed in section 
II, and is given by [20] 

S= (lnZ + ^\ -l(^!L\. XEtoMi (A8) 

where xeum is the topological susceptibility given by 
eq. (A4). One can explicitly see that the negative con- 
tribution is indeed present in the expression (A8) for the 
entropy. This term with the "wrong" sign in eq.(A8) 
is exactly proportional to the topological susceptibility 
(Al) in agreement with conjecture [14]. Furthermore, 
this term can be represented as a surface integral be- 
cause Q = 4__E entering (Al) is the topological charge 
density operator which is a total divergence. One should 
also emphasize that the entropy S as well as its surface 
term ~ xe&iM separately are gauge invariant observables. 
Also, while the term ~ Xeum can be represented as a sur- 
face integral, the entropy itself does not possesses such 
a surface representation. Finally, the entropy (A8) can 
be interpreted as the entanglement entropy because the 
only local observable is E, which is constant over space 
as shown in section II. It means that the measurements 
of E will be perfectly correlated on the opposite sides of 
the system [20]. 

The crucial observation for the present paper is as fol- 
lows. When the IR physics is properly treated, the en- 



tropy (A8) is obviously a positively defined function. Fur- 
thermore, as the theory under discussion is "empty" in 
the sense that it does not describe any physical propa- 
gating degrees of freedom in the bulk, one should expect 
that the entropy S must vanish in the infinite volume 
limit V — > oo. This expectation follows from the fact 
that the only dynamics in this system could be related 
to the so-called "edge states" which are localized at the 
boundary of the system but not in the bulk, similar to 
other topological field theories [7]. The only way this 
vanishing result could occur is the presence of a negative 
contribution which could cancel a conventional positive 
contribution present in (A8). In different words, the neg- 
ative contribution in (A8) is a must in order to produce 
an anticipated vanishing entropy in the infinite volume 
limit V — > oo. 

As we discussed previously, this contribution with a 
"wrong" sign can not be identified with any physical 
propagating degrees of freedom. Rather it is related 
to the tunnelling processes between different topological 
sectors as discussed in section II. These discussions again 
support our claim that the topological sectors must be 
included into consideration for self-consistency of the the- 
ory. Therefore, the additional terms they produce leading 
to the Topological Casimir Effect should be considered 
as physically observable quantities. As the last comment 
of this Appendix: though the term (A4) with a "wrong" 
sign is a gauge invariant contribution, its explicit compu- 
tation depends on a specific gauge-dependent technique 
being used. In particular, in the KS framework [18] this 
term is saturated exclusively by an unphysical ghost field, 
see explicit computations in [14]. Still, this term (A4) is 
physical as we argued above, and it can not be discarded 
on a sole basis that it is saturated by the artificial ghost. 

The main lesson of this Appendix is that there are 
extra contributions to the vacuum energy due to non- 
trivial topological features of the gauge fields, which do 
not have counterparts in scalar field theory. Therefore, 
the standard receipt (that the contribution to the energy 
and pressure due to the physical Maxwell photons could 
be obtained directly from expressions for massless scalar 
field by multiplying the corresponding scalar expressions 
by factor two) does not represent a complete description 
of the ground state in the presence of the gauge fields. 
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